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We analyzed 6 hours of data from the TAMA300 detector by matched filtering, searching for gravitational 
waves from inspiraling compact binaries. We incorporated a two-step hierarchical search strategy in matched 
filtering. We obtained an upper limit of 0.59/hour (C.L.=90%) on the event rate of inspirals of compact binaries 
with mass between 0.3Mq and WMq and with signal-to-noise ratio greater than 7.2. The distance of IAMq 
(0.5Mq) binaries which produce the signal-to-noise ratio 7.2 was estimated to be 6.2kpc (2.9kpc) when the 
position of the source on the sky and the inclination angle of the binaries were optimal. 



Introduction: Several laser interferometric gravita- 
tional wave detectors are now under construction. 
These include LIGO ffj, VIRGO §, GEO600 §, and 
TAMA300 j§. The TAMA300 detector has been devel- 
oped over the past five years. It is a power-recycled, 
Fabry-Perot-Michelson interferometer, which consists of 
mirrors that are suspended by vibration isolation sys- 
tems. The differential armlength changes in the 300m 
Fabry-Perot arm cavities caused by propagation of grav- 
itational waves are monitored by the Michelson interfer- 
ometer. 



The TAMA300 detector became ready to operate in 
the summer of 1999. Most of the designed systems (ex- 
cept power recycling) were installed by that time. The 
first test operation was done on August 6th 1999. The 
first long-term data were taken for three nights between 
September 17th and 20th 1999 j§ [§. The total data 
length amounted to about 30 hours, with the longest con- 
tinuous lock time of the interferometer lasting nearly 8 
hours. The strain equivalent noise spectrum is given in 
Fig. 1. The best sensitivity was about 3 x 10 _20 /v / Hz 
around 900Hz. The data of the first two days contained 



'deceased 



1 



many burst noises due to the instability of the interferom- 
eter system and it was very difficult to use those data for 
the gravitatioal wave event search. However, the quality 
of the data of the final day was good enough to perform 
the search for gravitational waves. 

Currently, TAMA300 is the world's largest detector in 
operation. This is the first time that long continuous data 
from one of the large-scale interferometer of the new gen- 
eration quoted above were taken. By analyzing the data, 
we expect to obtain useful knowledge of the property of 
long and continuous data. It is also important to develop 
tools for analyzing such data as a first step to treat the 
much larger amount of data which will be obtained in 
the near future. This is the first paper on data analysis 
of TAMA300, in which we report the first result of our 
search of gravitational waves from inspiraling compact 
binaries. 

Inspiraling binaries: Gravitational waves from inspi- 
raling compact binaries, consisting of neutron stars with 
mass ~ 1.4M Q or black holes, have been considered to 
be the most promising target for laser interferometers. 
These compact binaries can be produced as a conse- 
quence of normal stellar evolution of binaries. It has 
been also suggested that MACHOs || in our Galactic 
halo may be primordial black holes with mass ~ 0.5Af Q . 
If so, it is reasonable to expect that some of them are in 
binaries which coalesce due to the gravitational radiation 
reaction 0. 

Inspiraling compact binaries are nearly ideal sources 
for data analysis, since their waveforms can be theoreti- 
cally calculated to a very high accuracy. However, in this 
paper, we did not take account of the effect of spin an- 
gular momentum. Thus, we cannot deny the possibility 
that we lose signals from binaries with spin. g. 

Data Acquisition and Calibration: The main signals 
from the detector were derived from the feedback sig- 
nal that keeps the interferometer in resonance, which 
was digitized by an analog-to-digital converter in 16 bit 
depth with 20kHz sampling rate. The noise amplitude 
due to quantization by the ADC was estimated to be 
10~ 2 — 10~ 3 times smaller than the instrumental (elec- 
tronic) noise. A precise sampling clock was generated 
using a GPS-derived 10MHz reference signal. All of the 
data were recorded in the "Common Data Frame Format 
H" on tape archives. 

The strain of a gravitational wave h(t) and the volt- 
age signal V(t) are related in the frequency domain by 
h(f) — F(f)V(f), where F(f) is a response function. 
We measured the full open-loop response function in the 
observed frequency range before and after each continu- 
ous operation. During the observation, we continuously 
monitored drifts of the response function by adding a 
single frequency sinusoidal signal into the feedback sys- 
tem and observing its amplitude and phase at different 
points of the loop This method has an advantage 

that the observation is not interrupted by the calibration 
signal. The accuracy of strain measurement was evalu- 



ated to Ah/h ~ 1% within the observed frequency band. 
Using Monte-Carlo simulations, we also confirmed that 
systematic errors in the outputs of matched filter caused 
by calibration errors were negligible compared with noise 
fluctuations. The details of this method will be presented 
in a separate paper |p"l| . 

Matched filtering: We denote the strain equivalent one- 
sided power spectrum density of the noise by S n (f). In 
order to calculate the expected wave forms, which are 
called templates, we used restricted post-Newtonian wave 
forms of order 2.5, in which the phase evolution was cor- 
rectly taken into account up to the 2.5 post-Newtonian 
order, but the amplitude was evaluated by using the 
quadrupole formula. As for the 2.5 post-Newtonian phase 
evolution, we used formulas derived by Blanchet et al. 

When the gravitational wave passes through the inter- 
ferometer, it produces a relative difference AL between 
the two armlengths L. The gravitational wave strain am- 
plitude is defined by AL — Lh(t). The wave form h(t) is 
calculated by combining two independent modes of the 
gravitational wave and the antenna pattern of the inter- 
ferometer as 

hit) = A[h c (t — t c ) cos a + h s (t — t c ) sin a], (1) 

where t c is the coalescence time, and h c (t) and h s (t) are 
the two independent templates with the phase difference 
7r/2. To construct filters, we need the Fourier transforms 
of h c (t) and h s (t). They were computed directly by using 
the stationary phase approximation, the validity of which 
was established by Droz et al. The parameters to 

distinguish the wave forms are the amplitude A, the two 
masses mi, rri2, the coalescence time t c and the phase a. 
We did not include spins of the stars in the parameters. 

We denote the data from the detector as s(t). We de- 
fine a filtered unnormalized signal-to-noise ratio p after 
the maximization over a as 



p=V( 

(a,b) = 



s,h r ) 2 + (s,/i s ) 2 



df- 



Su(\f\) 



(2) 
(3) 



where a(f) denotes the Fourier transform of a(t) and the 
asterisk denotes the complex conjugation. This p has an 
expectation value v2 in the presence of only Gaussian 
noise. Thus, the signal-to-noise ratio, SNR, is given by 
SNR= p/V2. 

Analyzing the real data from TAMA300, we found that 
the noise contained a large amount of non-stationary and 
non-Gaussian noise whose statistical properties have not 
been understood well yet. In order to remove the in- 
fluence of such noise, we introduced a x 2 test of the 
time-frequency behavior of the signal |Q. We divide 
each template into n mutually independent pieces in the 
frequency domain, chosen so that the expected contri- 
bution to p from each frequency band is approximately 
equal. For two template polarizations h c (t) and h s (t), we 
calculate % 2 by summing the square of the deviation of 



2 



each value of p from the expected value fLa] . This quan- 
tity must satisfy the x 2 -statistics with 2n — 2 degrees of 
freedom, as long as the data consists of Gaussian noise 
plus chirp signals. However, there was a strong tendency 
that an event with large x 2 has a large value of p. Thus, 
by applying a threshold to the x 2 value, we can reduce 
the number of fake events without significantly losing the 
detectability of real events. For convenience, we renor- 
malized x 2 as X 2 /{% n ~ 2)- In the current analysis, we 
chose n = 16. This number was determined mainly by 
the limited memory of our computer. 

We searched the parameter space of 0.3M© < mi,rri2 
< 10M© with total mass less than IOMq. The low 
mass limit was chosen so that it covers the estimated 
mass of MACHOs as much as possible. In this parame- 
ter space, we prepared a mesh. The mesh points define 
the templates used for search. The spacing of the mesh 
points was determined so as not to lose more than 2 % of 
signal-to-noise due to the mismatch between actual mass 
parameters and those at mesh points. Using geomet- 
rical arguments, we introduced a new parameterization 
of masses that simplifies the algorithm to determine the 
mesh points ]T^] . 

The parameter space defined by using our new mass 
parameters turned out to contain 2057 templates in the 
present analysis. Although this number is not too large 
to perform matched filtering with a simple algorithm, 
when the noise power spectrum of data is improved, the 
necessary number of templates will increase. Thus, we 
introduced a two-step hierarchical search algorithm and 
developed tools of matched filtering which can be used 
to analyze future much larger data streams. 

two-step hierarchical search: The basic idea of the two- 
step search is simple. The templates for the first step 
search are chosen at the mesh points with coarse spac- 
ing. We perform the first step search with an appropriate 
threshold which must be chosen to be sufficiently low not 
to lose real events. When there are events that exceed the 
threshold, we search the points around them in detail. 

Before starting the two-step search, we evaluated the 
values of p and x 2 of all the data with a few selected 
templates to obtain an estimate of the background dis- 
tribution of p and x 2 , assuming all of these events are 
due to noise. Based on this distribution, we looked for 
several combinations of (p* , x* 2 ) so that there were only 
a few fake events which satisfy p > p* and x 2 < X* 2 si- 
multaneously. These thresholds became our preliminary 
second step thresholds. The first step threshold and the 
first step length of the mesh separation were determined 
to satisfy the following condition (CI): When p and x 2 
of an event at the second step mesh point satisfies a sec- 
ond step threshold, such an event must also be detected 
at, at least, one of the first step mesh points around the 
second step point with more than 98 % probability. By 
imposing this condition, it is guaranteed that the false 
dismissal rate of real detectable events due to the pres- 
ence of the first step search is kept to be small 16 . Note 
that the preliminary second step thresholds, (p* , X* 2 ), 



which were used to decide the first step thresholds, de- 
termined the detection probability of events with given 
amplitude. These preliminary second step threshold are 
given in Table 1(a). 

We performed simulations by adding signals to real 
data, and looked for combinations of the first step thresh- 
olds and the mesh separation which satisfy the criterion 
(CI). We chose one of them which minimizes the com- 
putation time in total. We introduced, at the first step, 
the x 2 threshold as well as the p threshold to select the 
candidates for the second step, which reduced the com- 
putation time for the second step. 

In order to reduce the computation time, we further in- 
troduced various techniques, some of which are explained 
in @. 

Here, we mention the cut off frequency f c of templates. 
The post-Newtonian wave forms must be terminated at 
some frequency roughly corresponding to the inner most 
stable circular orbit, at which the binaries are expected 
to transit to plunge orbit to coalesce. However, it is dif- 
ficult to determine, theoretically, the optimal value of f c . 
At the first step search, the choice of f c is important 
because the mesh separation is large. If the signal was 
located at a point far from the first step mesh points, the 
value of x 2 could become much larger than unity due to 
the difference of f c between the first step template and 
the signal. Thus, at the first step search, we searched 
for the value of f c which produced the largest value of p. 
On the other hand, at the second step search, we simply 
adopted f c = c 3 /(7 3/2 itGM) where M is the total mass, 
G is the Newton's constant, and c is the speed of light. 
This corresponds to the orbital radius of r ~ 7GM/c 2 . 

Results: Our analysis was done with 8 Compaq Alpha 
machines. Each of them can perform the FFT of data 
consisting of ~ 10 6 single precision real numbers with 
140MFlops. 

From the 30 hours of total data, we selected the data 
with better quality taken between 14:42 and 21:04 19th 
September 1999 (UTC time). The data consist of two 
separate stretches when the interferometer was contin- 
uously locked. An interval of 12.5 minutes around an 
unlocked part was not used for matched filtering. The 
total length of data used for matched filtering was 6.1663 
hours. 

In Figg we show a contour plot of the signal-to-noisc 
ratio as a function of the distance to the source and the 
total mass for equal mass binaries. We found that the 
signal-to-noise ratio becomes maximum around 1.6A/q 
binaries and decreases above this mass. This is because 
the noise level of the detector increases rapidly below 
several hundreds Hz. 

In performing the two-step search, the data were di- 
vided into data segments of 265.42 seconds with an over- 
lap of 29.9 seconds, which was longer than the longest 
template. The strain power spectrum density S n (f) was 
estimated using 530.84 seconds of data near the segment 
which was analyzed. 
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As a result of the second step search, we obtain p and 
X 2 as functions of masses and the coalescing time t c . In 
each small interval of coalescing time At c , we looked for 
an event which had the maximum p and which satisfies 
certain \ 2 threshold. In Fig.||, we show the numbers of 
these events which survived after the maximization as 
a function of SNR for various x 2 thresholds. The time 
interval At c was chosen so that the numbers of corre- 
lated events were small enough, but at the same time 
there was a sufficiently large number of events to investi- 
gate the properties of its distribution. In Fig.^J, we chose 
At c = 2 seconds. 

The computation time needed to perform the two-step 
search depends strongly on the first step threshold. If 
we had optimized the first step threshold by applying 
the above condition (CI), we would have obtained the 
threshold given in Table 1(b) in the case when the unit 
area of the parameter space of the first step search is 72 
times larger than that of the second step. In that case, 
the computation time would have been about 1.5 hours. 
However, we used a lower first step threshold as given 
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in Table 1(c) to obtain the result of Fig. 
vestigate the background distribution w 
statistical analysis. The computation time thus needed 
to obtain the results in Fig.|] was about 5.5 hours. In 
this analysis, we used a fixed set of the first step thresh- 
olds and spacing for all of the data. Thus, the number of 
events which satisfy the first step threshold, and there- 
fore total computation time, differs very much among the 
portion of data. 

Using this result, we estimate an upper limit on the 
event rate. In Fig|^, we show a fitting to the num- 
ber of events as a function of SNR 2 . The value of the 
X 2 threshold, 1.5, was chosen so that there was a 3.8% 
chance of rejecting real events in Gaussian noise. We 
chose the SNR threshold as SNR= 7.2. Thus, the ob- 
served number of events which exceed this threshold is 
2. Using this analytic fitting, we estimate the expected 
number of background events Nbg which is larger than 
SNR= 7.2 as Nbg = 2.5. Thus, using Bayesian statis- 
tics, and assuming uniform prior probability for the real 
event rate and the Poisson distribution of real and back- 
ground events, we estimate with 90% confidence that the 
expected number of real events which exceed SNR= 7.2 is 
smaller than 3.68 in this data M|. Thus, we obtained the 
upper limit of the event rate as 0.59/hour (C.L.= 90%). 
Note that we did not assume that the observed 2 events 
were real or not. (The largest SNR event for x 2 < 1-5 has 
SNR=9.0 and x 2 — 1-35. Although this event seems to 
show a small excess from the background distribution for 
X 2 < 1.5, it also seems to be absorbed in the background 
for x 2 < 2.5. ) 

From Fig. 2, we find the distance of a source that would 
produce SNR = 7.2 as 2.9kpc for mi = m 2 = 0.5M Q , 
and 6.2kpc for mi = m 2 = 1.4M Q respectively. In this 
paper, we estimate the distance in the case when the po- 
sition of the source on the sky and the inclination angle 
of the binaries were optimal. In order to evaluate the 



efficiency, we performed simulations. We added artifi- 
cial signals to the data and performed the same two-step 
search which derived Fig.|[ We then evaluated the detec- 
tion probability of artificial signals with the final thresh- 
olds determined in FigJ^. The result is given in Figjij. 
We found that TAMA300 could observe 1.4M events in 
several kpc. In other words, TAMA300 with the present 
sensitivity can give a direct observational limit of the co- 
alescing rate within several kpc. Here, it is important 
to note again that a new detector, TAMA300, produces 
data which are good enough to perform data analysis and 
to give upper limit of the event rate, even if it is sensitive 
only to events within several kpc. 

The sensitivity of TAMA300 detector is now being im- 
proved rapidly. With its goal sensitivity, 1.4M Q (0.5M©) 
binaries at distance D = 820kpc (D = 340kpc) will pro- 
duce a signal-to-noise ratio SNR=10. It is planned that 
much longer data with improved sensitivity and stability 
will be taken during the year 2000. By analyzing such 
data, we will be able to obtain a much more stringent 
upper limit to the event rate. And if we are very lucky, 
we may find a plausible candidate for a real gravitational 
wave signal. 
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FIG. 3. The number of events for x" < 2.5, X* < 1-5 and 
X 2 < 1-0 as a function of SNR 2 , and an analytic fitting to 
X 2 < 1.5. The fitting was determined between SNR 2 =35 
and 50. This fitting gives the number of events larger than 
SNR=7.2 and x 2 < 1-5 as 2.5. The X 2 < 2.5, x" < 1-5, and 
X~ < 1-0 correspond to 10~ 3 %, 3.8%, and 46% false dismissal 
rate in Gaussian noise, respectively. 
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FIG. 1. The strain equivalent noise spectrum of TAMA300 
on September 19th. The best sensitivity is about 
3 x 10 _20 /v / H^ around 900Hz. 
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FIG. 2. A contour plot of SNR as a function of distance 
and total mass of equal mass binaries in the case when the 
position of the source on the sky and the inclination angle of 
the binaries are optimal. 
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FIG. 4. The detection probability of artificial signals of 
1.4M0 binaries as a function of distance in the case when 
the final thresholds were (SNR*,x* 2 )= (7.2,1.5). 



TABLE I. (a) The preliminary second step thresholds used 
to determine the first step thresholds and the mesh spacing, 
(b) Optimized first step thresholds, (c) The first step thresh- 
olds used to derive the result of Fig. 3 and 4. In both cases, 
the first step mesh area is 72 times larger than that of the 
second step. 
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